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ABSTRACT 

Context. Asteroseismology is an efficient tool not only for testing stellar structure and evolutionary theory but also constraining the 
parameters of stars for which solar-like oscillations are detected, presently. As an important southern asteroseismic target, t Ceti, is a 
metal-poor star. The main features of the oscillations and some frequencies of t Ceti have been identified. Many scientists propose to 
comprehensively observe this star as part of the Stellar Observations Network Group. 

Aims. Our goal is to obtain the optimal model and reliable fundamental parameters for the metal-poor star r Ceti by combining all 
non-asteroseismic observations with these seismological data. 

Methods. Using the Yale stellar evolution code (YREC), a grid of stellar model candidates that fall within all the error boxes in the 
HR diagram have been constructed, and both the model frequencies and large- and small- frequency separations are calculated using 
the Guenther's stellar pulsation code. The^^ minimization is performed to identify the optimal modelling parameters that reproduce 
the observations within their errors. The frequency corrections of near-surface effects to the calculated frequencies using the empirical 
law, as proposed by Kjeldsen and coworkers, are applied to the models. 

Results. We derive optimal models, corresponding to masses of about 0.775 - 0.785 M and ages of about 8-10 Gyr. Furthermore, 
we find that the quantities derived from the non-asteroseismic observations (effective temperature and luminosity) acquired spectro- 
scopically are more accurate than those inferred from interferometry for r Ceti, because our optimal models are in the error boxes B 
and C, which are derived from spectroscopy results. 
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1. Introduction 

The solar five-minute oscillations have led to a wealth of infor- 
mation about the internal structure of the Sun. These results have 
stimulated various attempts to detect solar-like oscillations for a 
handful of solar-type stars. Solar-like oscillations have been con- 
firmed for several main-sequence, subgiant and red giant stars by 
the ground-based observations or by the CoRoT and the Kepler 
space missions, such as v Indi (Bedding et al. 2006; Carrier et al. 
2007), a Cen A (Bouchy & Carrier 2002; Bedding et al. 2004), 
a Cen B (Carrier & Bourban 2003a; Kjeldsen et al. 2005), n 
Arae (Bouchy et al. 2005), HD 49933 (Mosser et al. 2005),/? Vir 
(Martic et al. 2004a; Carrier et al. 2005a), Procyon A (Martic 
et al. 2004b; Eggenberger et al. 2004a; Arentoft et al. 2008; 
Bedding et al. 2010), t] Bootis (Kjeldsen et al. 2003; Carrier 
et al. 2005b), /3 Hyi (Bedding et al. 2001, 2007; Carrier et al. 
2001), 6 Eri (Carrier et al. 2003b), 70 Ophiuchi A (Carrier & 
Eggenberger 2006), e Oph (Ridder et al. 2006), CoRoT target 
HR7349 (Carrier et al. 2010), KIC 6603624, KIC 3656476 and 
KIC 1 1026764 (Chaplin et al. 2010), etc. Furthermore, the large 
and small frequency separations of p-modes can provide a good 
estimate of the mean density and age of the stars (UMch 1986, 
1988). On the basis of these asteroseismic data, numerous theo- 
retical analyses have been performed to determine precise global 
stellar parameters and test the various complicate physical ef- 
fects on the stellar structure and evolutionary theory (Thevenin 



et al. 2002; Eggenberger et al. 2004b, 2005; Kervella et al. 2004; 
Miglio & Montalban 2005; Provost et al. 2004, 2006; Tang et al. 
2008a, 2008b). 

t Ceti (HR 509, HD 10700) is a G8 V metal-poor star, be- 
longing to population II. Extensive analyse of this star have been 
performed by many scientists who have provided different non- 
seismic observational results (such as effective temperature T e ff 
and luminosity L), depending on the different methods used, i.e. 
interferometry and spectroscopy. Teixeira et al. (2009) detected 
solar-like oscillations on t Ceti, identified some possible existing 
frequencies, and obtained the large separation around Av = 169 
//Hz with HARPS. These seismological data will provide a con- 
straint on the fundamental parameters of r Ceti. Moreover, r Ceti 
will be one of the most promising southern asteroseismic targets 
of the seismology programme of Stellar Observations Network 
Group (Metcalfe et al. 2010). 

In this work, using a mixture of conventional and asteroseis- 
mic observed constraints, we try to determine modelling param- 
eters of r Ceti with YREC. The observational constraints avail- 
able to r Ceti are summarized in Sect. 2, while the details of 
the evolutionary models are presented in Sect. 3. The seismic 
analyses are carried out in Sect. 4. Finally, the discussion and 
conclusions are given in Sect. 5. 
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Table 1. Non-asteroseismic observational data of t Ceti. 



UDservauie 


value 


Source 


Effective temperature T e ff(K) 




i i *i 
(i) 




5525 ± 12 


(3) 


Luminosity L/L Q 


0.52 ± 0.03 


(4) 




0.50 ± 0.006 


(3) 




0.488 ± 0.010 


(2) 


Metallicity [Fe/H], 


-0.5 ± 0.03 


(1) 


Surface heavy-element 






abundance [Z/X] s 


0.0073 ± 0.0005 


(5) 


Radius R/R Q 


0.773 ± 0.024 


(5) 



References.— (1) Soubiran et al. (1998) , (2) Teixeira et al. (2009), (3) 
Pijpers et al. (2003a), (4) Pijpers (2003b), (5) this paper. 



2. Observational constraints 

2.1. Non-asteroseismic observational constraints 

The metallicity derived from observations is [Fe/H] =-0.5 ±0.03 
(Soubiran et al. 1998). The mass fraction of heavy-elements, 
Z, was derived assuming log[Z/X] « [Fe/H] + log[Z/X] , and 
[Z/X] Q = 0.0230 (Grevesse and Sauval, 1998), for the solar mix- 
ture. We can therefore deduce that [Z/X] s = 0.0068 - 0.0078. 
The radius, as an important parameter for constraining stellar 
models, was first measured by Pijpers et al. (2003a) using inter- 
ferometry. They determined the radius of t Ceti corresponding 
to 0.773 + 0.004(,'„ f .) ± 0.02(„,.)/? G . The measurement of the ra- 
dius was then improved by Di Folco et al. (2004) and Di Folco 
et al. (2007). Finally, Di Folco et al. (2007) determined the ra- 
dius R = 0.790 ± 0.005/? o . ln our work, we use a large value of 
radius R = 0.773 + O.O24/? which includes all the surrounding 
observational radius. 

The effective temperature and luminosity of t Ceti are both 
derived from spectroscopy (5264 ± 100K and O.52±O.O3L ), and 
by ensuring that we reproduce the measured radius ( 5525 ±12 
K, 0.500 ± O.OO6L ), using interferometry ( Soubiran et al. 1998; 
Pijpers et al. 2003a, 2003b). In addition the luminosity of a star 
can be obtained by combining our knowledge of the magnitude 
and distance. For t Ceti, the apparent magnitude V = 3.50±0.01, 
with the revised parallax, gives an absolute magnitude My = 
5.69 ± 0.01. Teixeira et al. (2009) derived a luminosity for t 
Ceti of L/L Q = 0.488 ± 0.010, using bolometric correction for t 
Ceti B.C. = -0.17 ± 0.02 (Casagrande et al. 2006) and adopting 
an absolute bolometric magnitude for the Sun of M ho i tQ = 4.74 
(Bessel et al. 1998). 

Using above different effective temperatures and luminosi- 
ties, we can obtain three error boxes, which error box A (5525 ± 
12 K, 0.50 ± 0.006 L ) are denoted by crosses, error box B 
(5264 ± 100 K, 0.52 ± 0.03 L Q ) denoted by triangles, and error 
box C (5264 ± 100 K, 0.488 ± 0.010 L Q ) denoted by diamonds, 
shown in Fig. 1(d), respectively. Meanwhile, we decided to in- 
crease all errors by a factor of 1.5, so that our calibration of the 
star is only weakly constrained by these values. 

All non-asteroseismic observational constraints are listed in 
Table 1. 



2.2. Asteroseismic constraints 

Solar-like oscillations of the G8V star t Ceti were detected by 
Teixeira et al. (2009) with the HARPS spectrograph. Thirty-one 
individual modes are identified (see Table 1 in Teixeira et al. 
2009). The large frequency separation is about Av = 169 jiHz. 



Table 2. Input parameters for model tracks. 



Variable 


Minimum 
Value 


Maximum 
Value 


6 


Mass M/M Q 


0.770 


0.795 


0.005 


Mixing length a 


0.8 


1.8 


0.2 


Initial heavy element 
abundance Z, 


0.001 


0.008 


0.0005 


Initial hydrogen 
abundance X, 


0.70 


0.75 


0.01 



Note. — The value S defines the increment between minimum and max- 
imum parameter values used to create the model array. 




3.76 3.74 3.72 3.70 3.745 3.732 3.718 3.705 
Log(Tj LogCU 

Fig. 1. (a) All evolutionary tracks in the HR diagram; (b) 
Evolutionary tracks falling in the error boxes from pre-main se- 
quence to main sequence; (c) Blow up the evolutionary tracks 
falling in the error boxes in the main sequence; (d) The selected 
models falling in the error boxes. Error box A (5525 ± 12 K, 
0.50±0.006 L Q ) is denoted by crosses, error box B (5264± 100 K, 
O.52±O.O3L ) denoted by triangles, and error box C (5264 ±100 
K, 0.488 ± 0.010 L ) denoted by diamonds, respectively. 

3. Stellar models 

We calculated many evolutionary tracks using Yale stellar evolu- 
tion code (YREC; Demarque et al. 2008) by inputting different 
parameters shown in Table 2. 

The mass range are M — 0.770 - 0.795 M with the incre- 
ment value 0.005 M . Initial heavy element abundance range are 
Z, (0.001 - 0.008) with the increment value 0.0005 and initial hy- 
drogen abundance X, (0.70-0.75) with the increment value 0.01 . 
Energy transfer by convection is treated according to the stan- 
dard mixing-length theory, and the boundaries of the convec- 
tion zones are determined by the Schwarzschild criterion (see 
Demarque et al. 2008 for details of the YREC). We set the mix- 
ing length parameter a=0.8-1.8 with the increment value 0.2. 
Using these parameter space, we created the model array. The 
initial zero-age main sequence (ZAMS) model used for r Ceti 
is created from pre-main-sequence evolution calculations. These 
models are calculated using the updated OPAL equation-of-state 
tables EOS2005 (Rogers and Nayfonov, 2002). We used OPAL 
high temperature opacities (Iglesias and Rogers 1996) supple- 
mented with low temperature opacities from Ferguson et al. 
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Fig. 2. (&)Xi values derived from Eq.(2), plotted as a function of age; (b)xl values derived from Eq.(3), plotted as a function of age; 
(c)xlc values derived from Eq.(6), plotted as a function of age; (d) |ro - 1| values plotted as a function of age. 



(2005). The NACRE nuclear reaction rates (Angulo et al. 1999) 
were used. The Krishna-Swamy Atmosphere T-t relation is used 
for solar-like star (Guenther and Demarque 2000). All models 
included gravitational settling of helium and heavy elements us- 
ing the formulation of Thoul et al. (1994). 



Figure 1(a) shows that many evolutionary tracks cover all 
possible evolutionary status of t Ceti. According to the above 
four error boxes, we select all the tracks crossing the error boxes 
shown in Fig. 1(b). We only choose to study main-sequence mod- 
els , which are shown in Fig. 1(c). Meanwhile, we use the mass 
and radius to estimate the large separation according to Eq. (1) 
( Kjeldsen & Bedding 1995; Miglio et al. 2009). Furthermore, 
using the temperature, luminosity, radius, and larger separation 
(refer to the values from Teixeira et al. 2009) as constrainst, we 
select the models of r Ceti provided in Fig. 1(d) as candidates. 



We now consider a function that describes the agreement be- 
tween the observations and the theoretical results 

5 Qtheo _ Qobs 
/= 1 i 

where C represents the quantities L/L G , T e jf, R/R G , and 
[Fe/H] s and large frequency separation Av, C' heo represents the 
theoretical values, and C bs represents the observational values 
listed in Table 1. The vector <rC° bs contain the errors in these 
observations, which are also given in Table 1 . We also decided 
to adopt a large error (all errors are increased by a factor of 1 .5), 
so that our calibration of the star is only weakly constrained 
by these values, which is not precisely determined. Figure 2(a) 



4 



Y. K. Tang et al.: Asteroseismology model of the metal-poor star t Ceti 



- 10 

i o 

o -10 
II 

f -20 

-30 

S -40 




3000 


4000 


5000 


6000 


10 





























M2 






-10 












-20 












-30 












-40 












3000 


4000 


5000 


6000 



Fig. 3. The difference between observed and best-fit model fre- 
quencies, according to the left term of Eq.(4). Squares are used 
for / = modes, diamonds for I =1 modes, triangles for / = 2 
modes, and circles for I — 3. Dotted lines show the power-law 
function, according to the right term of Eq.(4). 



presents the values x\ versus age t of selected models that are 
shown in Fig. 1(d). We find that we cannot select an optimal 
model from Fig. 2(a). From Fig. 2(a), we find that it is difficult 
to select an optimal model depending mainly on the non-seismic 
constraints and Av, which was estimated by simply scaling from 
solar value using Eq.(l). Hence, a detailed pulsation analysis are 
needed in the next step. 

4. Asteroseismic constraints of fundamental 
parameters 

Using Guenther's pulsation code (Guenther 1994), we calculate 
the adiabatic low-/ p-mode frequencies, the large- and small- fre- 
quency separations (Av n j = v nJ - v„_ u and 6v n j = v n j - v„_i,/ + 2, 
defined by Tassoul 1980) of all the selected models. We compare 
the theoretical frequencies with the corresponding observational 
frequencies using the function xl 



vf eo (n) - vf s (n) 



(3) 



where, N=31 is the total number of modes, and vp °{n) and 
v^ bs (n) are the theoretical and observed frequencies respectively, 
for each spherical degree I and the radial order n, where cr = 
2[iHz (Teixeira et al. 2009) represents the uncertainty in the ob- 
served frequencies and xl values, plotted as function of age, are 
shown in Fig. 2(b). 

Since existing stellar models fail to accurately represent the 
near-surface layers of the solar-like stars, where the turbulent 
convection take place, the systematic offset between the ob- 
served and model frequencies appears. Furthermore, this offset 
between observed and best model frequencies turns out to be 
closely fitted by a power law (Christensen-Dalsgaard & Gough 
1980; Kjeldsen et al. 2008; Metcalfe et al. 2009; Dogan et al. 
2009, 2010; Bedding et al. 2010; Christensen-Dalsgaard et al. 



2010). In other words, this offset increases with increasing fre- 
quency shown in Fig. 3. This power law can be expressed using 
the equation 



v„bs(n) - nv, heo (n) = ai[v ohs (ni)/v„ 



(4) 



where v b s are the observed frequencies of radial and non-radial 
order, v\, est = riv t i ieo (n) are the corresponding calculated frequen- 
cies of the best-fit model, and v max is a constant frequency cor- 
responding to the peak power in the spectrum, which is taken as 
4490 fiHz for t Ceti and r/, a/, and b a re parameters described in 
detail by Kjeldsen et al. (2008), (For a different spherical degree 
/, the values of r and a are denoted by r/ and a/, respectively.). For 
the Sun and a solar-like star, the exponent b = 4.90 is appropri- 
ate, as has been provn by many scientists. We use the Kjeldsen 
et al. (2008) prescription to correct the theoretical frequencies 
from near surface effects. 

According to Eq. (4), we can use the following equation to 
obtain the corrected frequencies of models: 

Vcorrectin) = riv theo (n) + ai[v ohs (n)/v max ] h . (5) 
We define the function^ in a similar way to Eq.(3) as 



I v c " rrec '{n) 



f. 



(6) 



The values of xl c , plotted as a function of age are shown in 
Fig. 2(c). From Fig. 2(c), we can see that the values of xl c are 
lower than xl and their lowest values correspond to model ages 
from 8 to 10 Gyr. We conclude that the optimal model corre- 
sponds to the lower values of xl and ro - 1. From Figs. 2(c) and 
2(d), we infer that only two models Ml and M2 can be accu- 
rately described by the observational constraints. The difference 
between the observed and uncorrected model frequencies of Ml 
and M2 are shown in Fig. 3. The uncorrected and corrected fre- 
quencies of the optimal models Ml and M2 and the observa- 
tional frequencies are shown in Table 3. 

To clearly compare all of the theoretical frequencies of the 
models with observational frequencies, we provide echelle dia- 
grams of models Ml and M2 in Fig. 4. An Echelle diagram is a 
useful tool for comparing stellar models with observations. This 
diagram presents the mode frequencies along the ordinate axis, 
and the same frequencies modulo the large separations in abscis- 
sae. From Figs. 4(a) and 4(d), it can be seen that the uncorrected 
theoretical frequencies are not closely in agreement with the ob- 
served frequencies. The corrected theoretical frequencies indi- 
cated by Eq. (5) fit perfectly the observation shown in Figs. 4(b) 
and 4(e). Because the observed frequencies of orders n are not 
consecutive and the values of v b s (n) are very close to those of 
v theo (n), we substitute the v, heo (n)/v max for v ohs (n)/v max . Hence 
Eq. (5) becomes 

Vcorrecin) = rv theo (n) + a[v theo (n)/v max ] b . (7) 

From Figs.4(b), 4(c), 4(e), and 4(f), it can be seen that corrected 
frequencies given by Eq. (5) and (7) respectively are uniform 
and reproduce the observed frequencies perfectly. Furthermore, 
we can use the function xl c to select the fitting model parameters. 
As we all know, the suitable model parameters correspond to the 
lowest values of xl c , which can be clearly seen in Fig. 5. From 
Fig. 5, we can conclude that the mass is in the range 0.775 - 
O.785M , a is in the range 1.6 - 1.8, Z, in 0.0065 - 0.0075, and 
Xi 0.73 - 0.75. Hence, the model parameters of t Ceti can be 
constrained to within these narrow ranges. Finally, we list the 



Y. K. Tang et al.: Asteroseismology model of the metal-poor star t Ceti 



5 



6500 
6000 

5500 

i — i 

x 5000 

^ 4500 

* 4000 

3500 
3000 



before correction 

i i i i | i i i i | i i i i | i i i i 

a=1.6 Z;=0.007 X,=0.740 M=0.775M, 



after correction by eq.(5) 



after correction by eq.(7) 



1=2 1=0 


1=3 I 


AC 


OV 


in 


O V 


Ml 


O V 


iAD 


ow 


in 


OtV 


An 


»? 


cm 


w 


All 


V 


ACi 


» V 


tc 


f 


ATJi 


» VT 


An i 


» ? 


ACi 


0» VT 


AD I 


01 V 


ACi 


O V 


A □■ 


V 


i An i i i 


, 0, vl 



-to 



N 
I 

i i 



6500 
6000 

5500 

5000 

4500 

4i 



3500 
3000 



1 II 1 1 1 

1 

o- i.d Z[-U.UU/ 


Mil 

1 1 

v — n ~i ac\ u 
A;-U. /4U M 


1 1 1 1 

-U./ /jM fl 


1=2 1=0 


1=3 I 


/ n v 

=i b) 

V / 


M 




■ 


IA 


T 




i 


T7 






» T 


■ 


Mi 


V 




j 






km 


• 




A 






A ■ 


§ ? 




□1 


e tv 




II 


■ Vt 




■ 


• 




n 






■ 

i i i i 1 i i i i 


l , , , , l 


1 1 1 1 



N 
I 

i i 



6500 
6000 

5500 

5000 

4500 

4' 



3500 
3000 



1 1 1 1 1 1 1 1 
1 

«— i c 7 — n nm 


i i i i i i i i 
i i 

V — Pi 7/10 U — 77^U 
\|-U./4U M-U.//jM 5 


i=2 1=0 


1=3 1=1 (c) 


H 


OV 


A] 


OV 


j| 


OV 


Mn 


T 


AD 


OV 


AD 


9 T 




H 


a n 


V 


k u 


• V 


A D 


T7 


£ ■ 


• 7 


a a 


e w ' 


a n 


a vt 


a ■ 


• v 


a n 


V 


A ■ 


V 


i i i i l\ in i i 


I i O i iVl i i i i 



50 100 150 200 
v modulo kv [/xHz] 



6500 
6000 

5500 

i — i 

x 5000 



before correction 



=1 



4500 

4000 

3500 
3000 



1 1 1 1 1 1 1 1 1 

a=1.6 Z ; =0.007 


| 1 1 1 1 | 

<:=0.750 U-- 


Mil 

0.785M E 


1=2 1=0 

AD 


1=3 I 

OV 


„(d) 


AD 


V 




Mi 


O V 




iAD 


ow 




4D 


Of 




AD 


»* 




m 


w 




AD ■ 


O V 




ifll 


» V 




AO 


» 




jDi 


» Vt 




AD I 


» » 




ACi 


Ot VT 




ACI 


a» v 




ACi 


O V 




ACi 


V 




i i i i 4 D i i i 


1 , 4), Vl 


i i i i 



N 
I 



6500 
6000 

5500 

5000 

4500 

4000 

3500 
3000 



50 100 150 200 
v modulo kv [jiiHz] 

after correction by eq.(5) 

i i i i i i i i i i i i "i i r | j | 



i i i i I i i i i I i 

0=1.6 Z=0.007 X,=0.750 M=0.785M, 



1=2 1=0 



14 
1 

AC 
J 

km 
i 

a ■ 
□i 
w 
j 
■ 
■ 



1=3 1=1 



• T 
• 

T7 
§ H 
TV 
» W 
I 



to 



N 
I 



6500 
6000 

5500 

5000 

4500 

4000 

3500 
3000 



50 100 150 200 
v modulo kv [jitHz] 

after correction by eq.(7) 

i i i i i i i i i i i i * i i r 



1 1 1 1 1 1 1 1 1 

0=1.6 ZpO.007 


I I 

X—0.750 U 


1 1 1 1 

=0.785M 6 


1=2 1=0 


1=3 1 


■.(()■ 


H 


OV 


AC 


OV 




m 


OV 






OT 




AD 


OT/ 




AD 


• T 




Ml 


O V 




A J 


O V 




ll 


• V 




A C 


T7 




A I 


• V 




A rj 


t> w 




a n 


• VT 




A ■ 


t V 




a n 


O V 




A ■ 


O V 




i i i i A iD i i 


1 , 0, ,vl 


i i i i 



50 100 150 200 50 100 150 200 50 100 150 200 
v modulo kv [jiiHz] v modulo kv [/xHz] v modulo kv \\iWi\ 

Fig. 4. Echelle diagrams for the optimal models Ml (upper panel) and M2 (lower panel). Left panel shows the case before applying 
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model parameters and characteristics of models Ml and M2 in 
Table 4. 



5. Discussion and conclusions 

Using the asteroseismic analysis and the empirical frequency 
correction for the near-surface offset presented by Kjeldsen et 
al. (2008) to correct our theoretical frequencies, we have derived 
the optimal model of t Ceti and now list our main conclusions: 
1. Using the latest asteroseismic observations, we have at- 
tempted to construct the optimal model of t Ceti. We have only 
considered the models Ml and M2 , which can closely describe 
the observations, as the optimal models. Furthermore, the model 
parameters of t Ceti have been constrained to within narrow in- 



tervals by the function xtc where the mass is in the range M = 
0.775 - 0.785M o , the mixing length parameter in the range a = 
1 .6 - 1 .8, the initial metallicity in the range Z, = 0.0065 - 0.0075, 
the initial hydrogen abundance in the range X, = 0.73 - 0.75, and 
the age in the range t = 8 - 10 Gyr. 

2. We have found that the results of the non-asteroseismic 
observations (effective temperature and luminosity) inferred 
from spectroscopy are more accurate than those derived from 
interferometry for t Ceti, because our optimal models are in the 
error boxes B and C derived from our spectroscopy results. 
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Table 3. The observational frequencies and the theoretical frequencies for model Ml & M2 before and after correction for near- 
surface offset, respectively. 
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Table 4. Final model-fitting results for r Ceti. 



Modelling parameters 


model Ml 


model M2 


Mass M/M 


0.775 


0.785 


Mixing length a 


1.6 


1.6 


z, 


0.007 


0.007 


x, 


0.740 


0.750 


Model characteristics 


Effective temperature T e ff(K) 


5409 


5387 


Luminosity L/L Q 


0.47985 


0.47612 


Log(g) 


4.53187 


4.53365 


Radius R/R Q 


0.78994 


0.79339 


(Z/X) s 


0.00753 


0.00749 


Age (Gyr) 


9.5 


9.5 


< Av > (iiHz) 


170.9222 


170.9106 


< Avi > (jiHz) 


170.8621 


170.8381 


< Av 2 > (±iHz) 


171.0555 


171.0332 


< Av 3 > (ptHz) 


171.5120 


171.4870 


< Sv Q2 > (jjHz) 


10.013 


10.111 


< 6v [3 > (jjHz) 


18.034 


18.136 


Model corrected paraments 


^0 


1.000302 


1.000264 


r\ 


0.9993002 


0.9993007 


T2 


0.9984142 


0.9984387 


r 3 


0.9984967 


0.9984996 


a Q 


-10.59438 


-10.32439 


a t 


-8.270579 


-8.092409 


a 2 


-6.517972 


-6.377440 


a 3 


-5.891401 


-5.639216 



